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The Casimir-Polder force is an attractive force between a polarizable atom and a conducting 
or dielectric boundary. Its original computation was in terms of the Lamb shift of the atomic 
ground state in an electromagnetic field (EMF) modified by boundary conditions along the wall 
^— ^ and assuming a stationary atom. We calculate the corrections to this force due to a moving atom, 

' demanding maximal preservation of entanglement generated by the moving atom-conducting wall 

^s^j , system. We do this by using non-perturbative path integral techniques which allow for coherent 

back-action and thus can treat non-Markovian processes. We recompute the atom-wall force for a 
conducting boundary by allowing the bare atom-EMF ground state to evolve (or self-dress) into the 
interacting ground state. We find a clear distinction between the cases of stationary and adiabatic 
motions. Our result for the retardation correction for adiabatic motion is up to twice as much as 
that computed for stationary atoms. We give physical interpretations of both the stationary and 
adiabatic atom-wall forces in terms of alteration of the virtual photon cloud surrounding the atom 
£^ , by the wall and the Doppler effect. 

> 

o 
o 
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The physical system studied in this paper is an atom in a polarizable ground state near a conducting wall. The 
interaction of the atom with the quantum electromagnetic field (EMF) vacuum, whose spatial modes are restricted 
by the wall with imposed boundary conditions, generates a force that pulls it toward the conducting wall (for general 
Oh, discussion see Ref. |lj). The details of such a force is important in any experiments and applications in which an 
■ atom is held near a surface by a trapping scheme using evanescent waves or magnetic fields. The atom- wall force is 
divisible into two parts. First, there is the electrostatic attraction that the atom feels toward its image on the other 
side of the wall, called the van der Waals (vdW) force. Second is a quantum mechanical modification of the vdW 
force first calculated by Casimir and Polder 0. They dubbed the quantum modification "retardation" of the vdW 
force, because its source is the non-instantaneous transverse EMF. Extensions of Casimir and Polder's results for a 
polarizable atom were later derived by many authors 3, 4], including for an atom in a cavity 51 and near a dielectric 
wall [H0- Closest in philosophy to what is done in this paper is the work of Milonni in Ref. @. There, the author 
computes the second order alteration of the EMF mode functions due to the presence of the atom, from which the 
ground state energy shift is the expectation value of the interaction Hamiltonian in the altered vacuum 24] . However, 
the author neglects time dependence in the mode functions and thus neglects effects due to Doppler shifts of the EMF 
modes. Recently, retardation correction of the vdW force has been demonstrated experimentally [fJEI- Verification 
of the Casimir-Polder force can be viewed as a demonstration of the entangled quantum behavior of the entire system, 
since it involves the dressing of the atom by the EMF vacuum. 

Although Casimir and Polder and others' calculations do treat the quantum entanglement in the system, analysis 
up to now has been restricted to stationary atoms. It has been assumed (wrongly, as we shall show) that such a 
method can also treat the adiabatic motion of the atom. Adiabatic motion means in this context that as the atom 
moves, it continuously shifts into the position dependent stationary dressed ground state on a timescale much shorter 
than the timescale of motion. Treatments assuming that the atom is stationary or is instantaneously static exclude 
correlations that are developed in the system during the motion. The key point is that the adiabatic and stationary 
dressed vacuum states are not the same. An example where this situation is encountered generically and dealt with 
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in depth is in cosmology, specifically, quantum field processes in an expanding universe For stationary systems 
a vacuum state is well denned at all times (due to the existence of a Killing vector), but not for arbitrary dynamics, 
especially fast motion. However, for slow dynamics, adiabatic vacuum states can be defined and renormalization 
procedures constructed [13jIi3jIi3- The adiabatic method we use here is similar in spirit (though not in substance, 
as our purpose is somewhat different from that in cosmology). To predict motional effects, entanglement in the 
evolution needs to be accounted for theoretically. We use the influence functional (IF) method here, which keeps 
track of full coherence in the evolution to derive the force between the atom and the wall while allowing the atom to 
move adiabatically. In the case of a stationary atom, our result is in exact agreement with the Casimir-Polder force. 
In the case of an adiabatically moving atom, we find a coherent retardation correction up to twice the stationary value, 
thus our coherent QED calculation will make verifiable predictions. This paper shows the derivation and explains 
the cause due to coherent back-action. Section II outlines the model and details of the calculation. The results for 
stationary and adiabatic motion are then given in Section III, and discussed in Section IV. 



II. MODEL AND APPROACH 



In contrast to obtaining the force via the gradient of the ground state energy shift, we obtain it through the 
expectation value of an atom's center of mass (COM) momentum. Our system consists of an atom placed near a 
conducting wall. We assume an initially factorized state of the atom in its ground state and the EMF in its vacuum. A 
path integral technique is used to derive the ground state- EMF vacuum transition amplitude of the evolving system. 
Inclusion of coherent back-action allows the system to self-dress 0, ^3 and preserves maximal entanglement in the 
non-Markovian evolution of an atom-EMF quantum system. The expectation value of the momentum operator is then 
computed. In the path integral, Grassmannian and bosonic coherent states are used to label the atomic and EMF 
degrees of freedom, respectively. The position and momentum basis are used for the atom's center of mass degree 
of freedom. The major approximation applied here is a second order vertex approximation. With the second order 
vertex, the propagator is partially resummed to all orders of the coupling constant. The result is a non-perturbative 
propagator which yields coherent long time dynamics |2Ct l2ll | . The mass of the atom and the size of its external 
wavepacket are kept finite throughout the calculation. Only at the end of the calculation do we allow the mass of 
the atom to go to infinity and its extension shrunk to a point, while retaining finite terms due to their effect on the 
dynamics. 

Highlights of the calculation are given in this section and details are given in the Appendices. In Section IIA 
the Hamiltonian and spatial mode functions that describe the system are introduced. In Section IIB the transition 
amplitude of the EMF vacuum with the atom in its ground state is calculated in a coherent state path integral, with 
an effective action expanded to second order in the coupling (equivalent to a second order vertex resummation), and 
semiclassically in the COM motion. The momentum expectation value and the retardation correction force is then 
calculated from the transition amplitude in Section IIC. 



A. The Hamiltonian 

The spinless non-relativistic QED Hamiltonian is given by [25| 

ff=i^ + J-(p-eA)* + eV(X) + H b . (1) 

The first term is the COM kinetic energy of an atom with mass M. The second term is the kinetic energy of the 
electron sitting in the transverse EMF. The third term is the potential energy of the electron around the atomic 
nucleus. The last term is the energy of the free EMF. After taking the dipole approximation, and restricting to two 
internal levels of the atom, the Hamiltonian in minimal coupling form becomes (see Appendix A of Ref. [22] without 
the rotating wave approximation) 

H = 217 + huJoS + S - + k J2 Wkb k b k + H n + H I2 = H + Hl (2) 

k 

The operators S± are the up and down operators of the atomic qubit and Uq is the atomic transition frequency. The 
operators bk and b^ are the EMF mode annihilation and creation operators, and Wk are the frequencies of the EMF 
modes. The two parts of the interaction Hamiltonian are 

H n = nJ2-^\PegS + +PgeS-}-[u k h k + ulhl} (3) 
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H n = hV] [u k ■ uib k bi + ul ■ ui(J k i + 2b[bi) + ■ ujb^bj]. (4) 

The vector p eg is the dipole transition matrix element, which is defined as p eg — (e\p\g) — — imo;o(e|r|g). The vectors 
Uk contain the photon polarization vectors ek and the spatial mode functions /k(X), i.e., Uk(X) = £k/k(X). The 

coupling constants are g = — -y/ 8 ^f c " and A = ^/ 47r ^ ac ', with a being the fine structure constant. 

In the presence of a conducting plane the spatial mode functions of the EMF which satisfy the imposed boundary 
conditions are the TE and TM polarization modes 

Uki(X) = yjk|, x Zsin(fc z Z)e ik "- x (5) 

Ti 



u k2 (X) = \j — -[fc||Zcos(fc z Z)-ik z k||sin(k z Z)]e Ik ii- x , (6) 



and their complex conjugates. 



B. The Transition Amplitude 

The transition amplitude between the initial and final coherent states with initial and final positions is given by 

i r t+T 

(X f ,{z kf },4> ri t + T\exp[-- J i/(s)ds]|X l ,{z k J,^;t). (7) 

The transition amplitude relevant to the atom- wall force is the amplitude that the atom moves from Xj to X/ 
without the emission of any physical photons. This is a very good assumption, since the probability for physical 
photon emission is extremely small |16| . The initial and final states are thus characterized by the atom being in its 
ground state and the EMF in vacuum, with arbitrary COM position states. The initial and final coherent state labels 
can be set to zero to reflect those states, although during the evolution the system evolves freely, and the motion of 
the COM is affected by recoil from emission and re-absorption of virtual photons, 

i f t+T 

Jf|X / ;t + r,X < ;t] = <X / ;t + r|exp[--y ff(s)ds]|X i; t). (8) 

Normally, a variational approach would be a sensible way to compute the functional integrals that make up the 
transition amplitude. However, since in this case both the anti-resonant as well as resonant rotating wave terms are 
included in the Hamiltonian (i.e., no RWA), the variational equations for the Grassmann variables will have bosonic 
sources even when the EMF is taken to be in the vacuum. We know from earlier work that when a Grassmann field 
variable has a bosonic source, the variational technique cannot unambiguously define the evolution of the Grassmann 
variable. A better way is to leave the transition amplitude as a discrete product of infinitesmal propagators. The 
necessary functional integrals can then be computed recursively. Details are in Appendix A. After the EMF and 
Grassmann path integrals are evaluated, the transition amplitude from the initial motional state X^ to the final 
motional state X/ (while keeping the same initial and final atomic ground state and EMF vacuum) is given to 0(e 2 ) 
vertex by 

/( rt+T [A/T-v-2 pS 2 

#Xexp i / + ip l drJ2~ e^ + ^ s -^u k (X(s)) ■ u*(X(r)) 



£-u k (X(s)).u k (X(s)) + 0(e 4 )lds) 



(9) 



where p 2 = (g|p 2 |g) is the ground state expectation value of p 2 . 

A semi-classical approximation to the transition amplitude Eq. 10 is obtained by evaluating the action along its 
classical path. This will neglect the fluctuation terms of order O(jj). The classical path is the straight line path plus 

terms of order O(fy), 

X c (s) = X, + ^^(s -t)+o(^)= X°(s) + O(^) . (10) 
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Evaluating the transition amplitude along that path gives 
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/ m \ 3 ^ 2 ( r t+T rMX° 2 r s p- 2 

K[X f ;t + r,X t -t]=^^-j exp |i jf + ip 2 jf dr £ JL e -(^+^)(- r ) Uk(X »( s )) . u k (X°(r)) 



2 -<(X°(s)) ■ u k (X°(s)) + 0(e 4 /M) 



ds 



en; 



Using the spatial mode functions of Eqs. (|5I6|I in the above gives the semi-classical transition amplitude in the presence 
of a conducting wall (see Eq. iJHIJ)- 



C. momentum expectation and force 



Given the above expression for the transition amplitude and an initial center of mass wavefunction for the atom, 
^(P), the momentum expectation and the force on the atom (the time derivative of the expectation momentum) can 
be computed. The momentum expectation is 

,Aw x ft f dPf 

™ + ^ = nJ Wf 

with the normalization factor 



,3 P / 



dXidX^ K[P /; t + r|X i; t] *(X 4 )**(X0 X*[P /; t + r|X^t], 



(12) 



N 



J J dXdX^ X[P /; t + r|X;t] *(Xi)**pCj) ^[P /; t+r|X^t]. 



(13) 



The initial wavefunction can be taken to be a Gaussian centered at (R, Pq) with the standard deviations (cr, 1/ct). 
Such a choice will allow for the possibility that the atom and the wall are moving toward or away from one another. 
Following the line of calculation detailed in Appendix B, a momentum moment generating function is computed in 
the limits M — > oo and a — > such that jj — > V and a 2 M — > oo (see Eq. (IB12I) ). From the generating function the 
momentum expectation value can be computed, 



P(t + r) 



h dZ{3) 



iZ(0) dJ 



(14) 



j=o 



In the above limits 
P(t + r) = P 



2iA 2 h ^ k z cos 2 fl f t+T dg e _ 2ikz . (R+v(s -t)) 



L 3 ^ 

k 



g 2 Pz^ 

L 3 



E 



k z cos 2 1 



t+T 



<!s ^ dr e - ikz -( 2R + v ( s+1 '- 2t )) 



-i(oj k +wo)(s-r) _ i(w k +wo)(s— r) 



(15) 



The momentum depends on the position and velocity only through the distance from the wall and the velocity toward 
or away from the wall, so motions parallel to the wall have no effects. Define R = e z ■ R and v = e z • V, with e z 
defined as positive away from the wall. Taking the time derivative of the momentum expectation value will give the 
force that is exerted on the atom by the transverse EMF in the presence of the wall. Doing so, as well as applying 
the Thomas-Reiche-Kuhn sum rule, 



A 



2 2 
gPz 

UJ 



(16) 



and rewriting in terms of the static ground state polarizability, a Q , the force is 



F c (R,v,t + r) 



27ria bwo k z cos 2 6 



L 3 



-2ik z -(R+Vr) 



k 



iraohujQ x-^ k* cos 



L 3 



k 



2 L 



t + T 



ds 



-ik z -(2R+V(r+s-t)) 



-i(w k +w )(t+r-s) _ i(w k +u )(t+r-s) 



.(17) 



The subscript "c" is a reminder that the force calculated from the transverse field is the retardation correction to 
the electrostatic force. Inspection of the force reveals that it is a sum over recoil momenta weighted by amplitudes 
which depend on the distance of the atom from the wall and the velocity of the atom. As will be discussed in Section 
IV, the recoil momenta come from virtual photon emission and re-absorption. In that sense the net force reflects an 
interference phenomenon, since it is the net sum of many different possible virtual processes. 
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FIG. 1: This plot shows the value of the atom-wall force at R = 3000 vs time in atomic units. The spike at r = 6000 is the 
time at which a photon emitted at r = will have just returned. Before r = 6000 the force is experiencing transient behavior, 
and afterward it rings down to the stationary atom value. 
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FIG. 2: This plot shows a snapshot of the coefficient of the behavior of the the atom- wall force at a time r = 6000 in 
atomic units. The location of the spike at R = 3000 corresponds to the location at which a photon emitted at r = will have 
just returned to R = 3000. At locations R < 3000 the force has begun to asymptote to its steady state behavior, and those 
at R > 3000 are still experiencing transient behavior. The inset image is a magnification near the wall. The dotted line is the 
coefficient of the i dependence of a stationary atom. 



III. RESULTS 



A. Stationary Atom 

If the atom is stationary, then setting v = gives the retardation force to be 
F(°)(R,v = 0,t + r)= - ^^l^e-« 



L 3 '-7? 

k 



3 \r ™ c 2 a r t+r 

(18) 



Trapfajg K cos 2 8 f +T ^ e _ 2ikz . R 



c -i(w k +wo)(t+T-s) _ c i(w k +wo)(t+T-s) 



k 

Combining the correction force with the electrostatic force gives the total force on a stationary atom, 

F sa (R, t + r) = -e z + (R, v = 0, t + r) . (19) 

The stationary atom force exhibits a transient behavior when the atom first "sees" itself in the wall. Then, on a 
timescale of several atom-wall round trip light travel times it asymptotes to a constant steady state value. The 
transient behavior is plotted in Fig. (1) and Fig. (2) for an optical transition frequency in an alkali atom. The steady 
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state value of the stationary atom-wall force can also be determined analytically to be 

P sa (R,r»2R/c) = -e 2 ^0-e^fAy f dfc (20) 
V ' 1 ' ~ 8R 4 4tt \dRJ J kc + ujQ 2kR ' V ; 



which can be simplified to 



'^>>™-^{k)'*r <2i) 

From Eq. I|21[l the potential which a stationary atom feels is easily found to be 

V~m = -^(£) a i f^e-^, (22) 



with asymptotic limits 



8tt VdR/ R 



U sa (R)^-^^ forR«^ 
U sa (R) -> -^#- c ^ forR»^ 



(23) 



which exactly reproduces the results of energy gradient approaches. Although the results are the same as those 
previously derived, the interpretation behind how the results are obtained is different. The energy gradient approach 
can be described as a kinematic approach since the atom-EMF system is assumed to be held static in its entangled 
dressed ground state. The self-dressing approach used here, on the other hand, allows the atom-EMF entanglement 
to evolve dynamically. That is, the atom and EMF system, beginning in a factorized state, evolves into a stationary 
dressed state (i.e. it self-dresses). When the atom is stationary the two forces match because after some time to 'get 
acquainted', the self-dressing atom does indeed evolve into the stationary dressed state. It should be stressed that the 
agreement between the results of the two methods demonstrates the coherence of the self-dressing method as applied 
here. 



B. Adiabatic Motion 



We now show that the prediction given by the self-dressing method of the retardation correction force for a slowly 
moving atom differs from the energy gradient prediction |26| . The key difference is that as a moving atom and the 
EMF get acquainted, they evolve into an entangled dressed state which is different from the stationary atom dressed 
state. The reason for the difference is the Doppler shift of the EMF modes in combination with the presence of the 
wall. We will discuss this point in more detail in Section IV. 



1. adiabatic evaluation 



The retardation force for a moving atom can be determined from Eq. Ijl7|) by applying a separation of short time 
scale dynamics from long time scale dynamics, and how they affect each other. The adiabaticity condition is applied 
here in the same way that it is applied in standard methods for determining the dipole force on an atom in a laser 
beam j23|. There, assuming that the atom's position is constant on short timescales, the optical Bloch equations 
are solved for the steady state values of the internal state density matrix elements. On long time-scales the matrix 
elements are replaced by their steady state values and put into the Heisenberg equation of motion for the atomic 
COM momentum. Such a procedure is justified when the internal and external dynamics evolve on vastly different 
timescales. The analogous separation here will be of the short timescale describing the self-dressing of the atom-EMF 
system and the long timescale describing the motion of the atom. 

In order to be explicit about the timescale separation let us first rewrite Eq. I|17[) with the definition x = s — t, and 
remember that t is the time at which the atom-EMF system begins to evolve from a factorized state, 



F c (R t ,v t ,r) 



27ria 5a;g ^ kz cos 6 p -2ik z -(R t +v t r) 
L 3 l-j Wk 

k 



k 

iraghajQ v k 2 cos 2 i 



L 3 



£ 



Jo 



j x e -ik z -(2R t +V t (T+x)) 



-i(wk+u>o)(T— x) _ i(w k +U[))(T-x) 



(24) 
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so that the short timescale dynamics (parameterized by r and x) is explicitly separated from the long timescale 
dynamics (parameterized by t) on which R t and V t evolve. An adiabatic evaluation of the retardation correction for 
a moving atom can be extracted from a Taylor series expansion of Eq. Q24J1 , 



F c (R t ,v t ,r) = V^F(")(R t ,v t = 0,r) (25) 

n=0 



where n denotes the nth derivative with respect to velocity. The Taylor series expansion is an equivalent representation 
of the LHS as long as the RHS converges. Each function F( ra '(R t ,v t = 0, r) exhibits a transient behavior while the 
atom first "sees" itself in the wall (during times r ~ — ) and asymptotes to steady state behavior on a timescale of 
several round trip light travel times. The adiabatic approximation is applied at this point by replacing each function 
F( ra )(R t ,v t = 0,t) by its asymptotic behavior 

F(«)(R t) v t = ! r)^F(»)(R t) r) = _,/_l^(A_Y +3 f sin(2 fc R t ) 
v ' ; ss v z 2 n Air \dR t J J kc + LU 2kR t 

which means replacing the Taylor expansion, Eq. I|25|) . by its steady state form, 

oo n 

F c (R t ,v t ,r)^Ff(R t ,v t ,r) = ^^FW(R t ,r) (27) 

This step is analogous to replacing the internal state density matrix by its steady state value in adiabatic computations 
of the dipole force on an atom in a laser beam. Replacing the Taylor expansion by its steady state behavior is adiabatic 
because it assumes that the expansion terms asymptote to their dressed state form on a timescale much shorter than 
the timescale on which either the position or velocity of the atom changes. More specifically, for the change in position, 
the adiabatic condition means that during a round trip light travel time the atom- wall distance has very little relative 
change, v 2 ^ << R, which is equivalent to the condition that the atomic velocity be non-relativistic, 

I « \ . (28) 

Similarly, the adiabatic condition for the change in velocity is that it has very little relative change during a round 
trip light travel time, ^f 1 ^ << V, which can be restated as the net force not changing the kinetic energy of the 
atom much during a light travel time, 

{F net V)~ « ^A/V 2 , (29) 

c 2 

since F ne ± v is the power that the net force puts into the atoms mechanical motion. Both conditions are satisfied in 
typical experimental setups. 

Note that rather than tending to a constant steady state value, the terms in the Taylor expansion, Eq. I|26[l. 
asymptote to steady state polynomial time dependence, the source of the polynomial time dependence being the 
k z • V Doppler shift term in the exponents of Eq. 124|) . In distinction to the stationary atom case those polynomial 
time dependencies will lead to non-zero partial time derivatives as well as the convective changes due simply to motion 
the of the atom 

d fdR d dv d d\- , on , 

d~s Fc - (, ds"dR + ds~dv + d~s) Fc - (30) 

The differential change in F c can then be split into two parts, one coming from the convective change and the other 
from the partial time derivative, 

dF 

dF c = ^ c \ convectme + ds^f. (31) 

The convective differential change is the differential change in the force not including any short timescale time 
dependence, in other words, the steady state expression at r = 0, 



dF r = dF^. 



d, f: (32, 
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with, from Eq. Q27[l. dF^ s | T _ = dFc * 1 . The behavior of the force on long time-scales is computed by integrating the 
differential change from an initial time at which v = up to the final time, 



d Sl ^( Sl ), 



(33) 



where it has been substituted that P c (to) = Fc (to) (since v = at to). A similar analysis for the differential of the 
first partial time derivative gives, 



if— 

V 9s 



dF s c l 
~ds~ 



T = 



from which, 



ds 



(si) 



2 dR* c (Slj 



ds \ ds 



Sl d 2 F c 
ds 2 -^r(s 2 ). 



9s 2 



(34) 



(35) 



Carrying on similar analysis (and rewriting in terms of the zeroth order expansion term) leads to the general expression 



d n F c 
ds" 



(Sn) 



v Sl v S2 ..v s „ d n 
2" dR" 



gn+i F 



F(°)(s„)+/ ds„ +1 -^_ TI £( Sn+1 ) 



(36) 



Concatenating Eq. I|33|) with Eqs. (|36[) leads to an expression for the retardation correction force which is the sum of 
a series of imbedded integrals, 



p c (t)=Pi°)(t)+jf d S1 ^A F (o) (Sl)+ jr dsl ^ si dS2 ^a_ F (o) (S2 



(37) 



This result could have been written down directly since it has a straightforward interpretation of being the sum of 
the integrated effects of each of the partial time derivatives. Each term in Eq. I|37[l can be evaluated by making a 
change of variables from time to position with the identity v = dR/dt. For example, the first term gives, 



F(°>(R t )-Fi°)(Ro) 



and further terms give, 



dsi / ds 2 



ds, 



v Sl v S2 ..v Sn d T 



-F(°)(s 



2™ dR™ c v "'" 2 r ' 
Substituting these into Eq. (|37|l gives a geometric series with the result 



F(°)(R t )-F(°)(R ) 



(38) 



(39) 



f c (r) = f^(r) + j2(1) (f (o) (r)-f (o) (r : 

n=l ^ / \ 

= 2F(°)(R)-F(°)(Ro), 



(40) 



where R Q = R(to) is the distance from the conducting wall at which the atom was originally at rest. The force F^ ) (R) 
is the stationary atom retardation correction to the vdW force. 



2. force and potential 

Inspection of Eq. I|4U|) shows that if the atom is released but remains stationary, then the retardation force will 
be the stationary atom value. On the other hand if the atom is released infinitely far from the conducting wall and 
moves in toward the wall, then the retardation force near the wall will be twice the stationary value. At a finite 
initial distance the retardation force will vary between these values. The force in all cases will depend only on the 
position. Thus the atom still moves as if it were in a conservative potential and the potential it feels depends on 
where it started. 
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Combining the retardation correction force with the electrostatic force and simplifying as in Eq. i|21[l gives the 
atom-wall force to be 



.(R) 



a hujQ ( d \ 1 



dR ) R J x 



dx 



-2Rx/c - "ohwg ( d 



4tt 



dr 



dk sin(2kr) 



kc + u>o 2kr 



(41) 



The first term is the stationary atom-wall force and the second term is a residual force which pulls the atom back to 
its original point of release. The force can easily be turned into the potential which the atom feels: 



t^am(R) 



a hujQ ( d \ 1 



8tt \dR J R Jq x 



dx 



-2Rx/c 



* hug / d \ 2 f°° dk sin(2kr) 
4-7T \ dr / J kc + luq 2kr 



(42) 



Since the first term in the potential is the stationary atom-wall potential, in the regions near and far from the wall it 
will have the expected inverse powers of distance dependence, as shown in Eq. 123(1 . The second term is the residual 
potential due to the motion. 



IV. DISCUSSION 



A. Physical Interpretation 

In the energy gradient approach, one interprets the force between a polarizable atom and a wall as arising from the 
Lamb shift in the atomic ground state energy. Spatial variation of the ground state energy is expected to generate a 
force which pushes the atom to lower energy positions, but the mechanism for such a force is not given explicitly. In 
the final analysis, since the only players in the full system are the atom and the EMF field, such a force must come 
from the emission and reabsorption of photons. Our approach provides an interpretation of how a net force arises 
from the emission-reabsorption processes in the presence of a boundary. 

The connection between the Lamb shift calculation and our calculation is the dressed ground state of the atom, 
which is the true ground state of the full Hamiltonian. Expanded in the free (or bare) Hamiltonian basis, the dressed 
ground state is a quantum superposition of bare atom-EMF states, and is often described as an atom surrounded 
by a cloud of virtual photons which it continually emits and reabsorbs. In the energy gradient approach, the atom- 
EMF is assumed to always be in the stationary dressed ground state. By contrast, in our approach a bare state is 
allowed to evolve quantum mechanically into the dressed ground state. The difference between these two is crucial to 
understanding how the coherent QED correction comes about. By allowing the atom-EMF to evolve into a dressed 
ground state we leave open the possibility that the motion of the atom can affect how closely to the stationary dressed 
ground state the system evolves. Or in the language of the virtual photon cloud, the distribution of virtually occupied 
modes is allowed to differ from the stationary atom case. 



1. stationary atom 

Even without motion, the atom's virtual photon cloud is altered by the presence of the wall. For a perfectly 
conducting wall, the TE and TM spatial mode functions of the EMF are given by Eqs. JSJ-EJ- Those mode functions 
are determined by solving the wave equations with the given boundary conditions on the wall, and are constructed 
by linear combinations of plane wave modes. The creation and annihilation operators of the TE and TM EMF modes 
(bT, b) are thus combinations of the creation and annihilation operators of plane wave modes (a^, a) moving toward and 
away from the wall. Inspection of the Hamiltonian and the propagator shows that it is emission followed by absorption, 
which is the source of the force. In the interest of finding a physical interpretation, one can think of virtual processes 
in the presence of the wall in terms of plane waves. Then the emission-reabsorption of a wall-constrained mode is: 

b k bt Uk (X) ~ (a k e ikX -a_ k e- ikX )(a k e- ikX -aL k e ikX ) 

- a k a[ + a_ k al k - a_ k a k e~ 2ikX - a k a! k e 2ik ' x 

The first two terms are emission-reabsorption of the same photon and contribute no net momenta to the atom. The 
second two terms are emission of one photon and reabsorption of the reflected photon. Each of those contributes a 
2k 2 momentum to the atom. The effect of those processes on the force can be seen explicitly in Eq. I|24|) . The first 
term in Eq. (|24|l originates from the iff 2 interaction and the second terms from the Hn interaction. In both terms, 
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the sum over wavevectors is a sum over emission followed by reflected absorption processes, with each contributing a 
2k z momentum. Thus, the presence of the wall alters the atoms virtual photon cloud by reflecting some of the modes. 
The process of emission and reabsorption puts the photon cloud into a steady state distribution with the net effect 
on the atom of a retardation force. 



2. moving atom 

Once the stationary retardation force is understood in terms of the wall effect on the virtual photon cloud, the 
modification of it for an adiabatically moving atom can be interpreted as part of the Doppler effect. The effect is 
easiest to explain in the reference frame of the atom, in which it is the wall which will be moving toward or away 
from the atom. Then, as in the stationary case, the virtual photon cloud will be altered by reflection off the wall. 
However, in the case of the moving wall, the reflected photons will be Doppler shifted due to the walls motion. In 
the language of the virtual photon cloud, the distribution of photons around a moving atom will be Doppler shifted. 
This shift builds up in the photon cloud much like charge in a capacitor connected to a loop of wire in a changing 
magnetic field, and it can only be discharged through absorption into the atom. The net effect, over the retardation 
force, will be to push the atom against such built up Doppler shift, back to its original point of release. 

B. Prospects For Experimental Observation 

1. Reflection From An Evanescent Laser 

A situation in which the motional modification of the retardation correction will be important is for the reflection 
of cold atoms off the evanescent field of an otherwise totally internally reflected laser beam. For example, in a recent 
experiment by Landragin et. al. flp| . cold alkali atoms are dropped onto a crystal with an evanescent wave running 
along the surface. The atom-wall interaction pulls the atoms towards the wall. The dipole potential of the evanescent 
wave, on the other hand, causes a repulsion of the atoms from the crystal. The combination of those two creates 
a barrier through which some fraction of the atoms tunnel and the rest reflect back out. The authors measure the 
fraction of reflected atoms versus the barrier height. As the barrier height is lowered it will at some point drop below 
the energy of the incoming atoms. At that point, all the atoms will be able to classically roll over the barrier, and no 
atoms will be reflected. The evanescent laser power required to reach that barrier height depends sensitively on the the 
atom-wall attraction. By comparison of measurement with theory, the authors show that the electrostatic attraction 
alone does not accurately predict the threshold laser power. They show that the prediction of a retardation corrected 
force is closer to the measured value. When we combine the motional modification to the retardation correction we 
are able to make a further modified prediction for the threshold. The calculations done in this paper are for a perfect 
conductor, not a dielectric boundary, so the modifications predicted here should not be applied directly to the case of 
a dielectric boundary. However, a general statement can be made that a coherent QED correction will cause a lowered 
prediction for the threshold laser power, since it will tend to decrease the atom-wall attraction. If one naively applies 
a dielectric factor to our result for the conducting plate to compensate for the difference, the present prediction for 
the threshold energy (14.8 V) is close to the measured value (14.9 T). Extension of the present work to a dielectric 
wall is ongoing. 

2. Transmission Between Parallel Plates 

Another experiment which has been able to observe the retardation of the van der Waals force involves a stream of 
ground state atoms passing between two plates 9J. Due to the attraction of the atoms toward the plates, some of the 
atoms fall onto and stick to the plates. The fraction of atoms that pass through the gap depends on the atom-wall 
potential. By measuring the opacity (fraction of atoms that do not pass through) for different gap widths, the authors 
probe the attractive atom-wall potential. This experiment holds less promise of observing a coherent QED correction 
to the retardation, than the previous example. The reason being that in this experiment the atoms first come into 
interaction with the walls at a distance of only a few resonant atomic wavelengths. The atom and EMF thus do not 
have as much motion over which to develop a coherent effect. Within that caveat, a general prediction can be made 
that the coherent correction will tend to decrease the opacity. 
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C. Conclusion 

Our result exactly reproduces the Lamb shift result for a stationary atom. For an adiabatically slowly moving 
atom, a correction due to the Doppler shift is found. Our result for the retardation correction for adiabatic motion 
is up to twice as much as that computed for stationary atoms. Agreement with the energy gradient result in the 
stationary atom case shows that our non-perturbative approach captures the effects of entanglement which we sought. 
The physical interpretation is that the atom-EMF system evolves from an initially factorizable bare state into the 
interacting Hamiltonian ground state, which is an entangled state in the free Hamiltonian basis. This process is 
known as self-dressing. The correction for a slowly moving atom shows how our approach can go beyond Lamb shift 
calculations. The correction is due to the Doppler shift in that the virtual photon cloud which dresses the atom is 
shifted. 

APPENDIX A: RECURSIVE CALCULATION OF EFFECTIVE ACTION 

The Hamiltonian is given in Eq. J5J). The evaluation of the transition amplitude as a path integral begins with 
slicing it into infinitesimal steps. A single infinitesimal step transition amplitude for initial EMF vacuum and atomic 
ground state (i.e. the initial EMF and Grassmannian labels set to zero) is, 



(Xi, {z lk }, Vi; t + e| exp[--He]\X , {0 k }, 0; t) 

n 

"iM(Xi -X ) 2 e 
= 6XP [ 2^ 

'iM(X x - X ) 2 e 
= 6XP 2^h 



E7 _ gie f v-^ - - ^ e t t 



kl 



Al + X! V'l^lk-Bike + ^2 Zlk^llC'ikl 



k< 



kl 



(Al) 
(A2) 



With the obvious definitions of A\, -Bike, and Ci jk i. The first infinitesimal step transition amplitude, ea i|All) . can be 
used to derive the 2 infinitesimal step amplitude: 



(X 2 , {z k2 }, ^ 2 ;t + 2e| exp[-2-Fe]|X Q , {0 k }, 0; t) = 

d^(Xi)d^(zi)d^(V>i)(X 2 , {z 2 k}, V>25 1 + 2e| exp[-^-JJe]|Xi, {zi k }, ipi; t + e) 

n 



x (Xi,{2i k },Vi;t + e|exp[--ffe]|Xo,{O k },0;t) (A3) 



The result is: 



(X 2 , {z 2k }, t + 2e| exp[-2-He] |X 0j {0 k }, 0; t) 



d/x(Xi) exp 



A 2 + 2^'>p2Z2kB 2ke + 2^ z 2k2 ; 21 C 2kl + 2^ - 



kl 



2e 2 h 
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For definitions of the coefficients see Eq. I|A6(1 with n — 2. The 2-step transition amplitude can be generalized to an 
n-step transition amplitude: 



it 

(X„, {z„k}, i> n \ t + ne\ exp[— — 



/n r n 

k-B nke + ^ Cnkl + / ] 
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(A5) 



with the finite difference equations: 



A n = A n -y - ie — (u^ k • u nk ) - ie Y] -^L(p g0 • u nk )B n _i. ko + 0(e 2 ) 



CJ k 



ko 



(A6) 
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In the continuous limit those become first order differential equations with the following integral solutions: 
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The transition amplitude of Eq. l|A5Jl can be further simplified by setting the final EMF and atomic states to vacuum 
and ground, respectively. The transition amplitude is then: 



t+T 
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ds 
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The equations for B(s) and C(s), Eq. (|A10|) and Eq. (|A11|I . are Volterra type integral equations. Their solutions are 
infinite Born series in orders of the coupling. Approximations in the above coefficients are approximations in the basic 
vertex. To (9(g 2 ): 



A(t+r) = -i / T ds]T-[u k (X(s)).Uk(X(s))]- / ' ds / dr£ — e- i ^+-)( s - r )[u k (X(s)).p go ][u k (X(r)).p og ] 
Jt k J t Jt ^ ^k 

The transition amplitude with an 0(g 2 ) vertex is thus: 
(X /; t+r|exp[-^y + ffds]|X i; t) = 



DXcxp I i 
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In the above transition amplitude the polarization mode functions are dotted with the dipole vector of the atom. 
The direction that the atom's dipole vector takes will depend on the quantization direction chosen for the atom's 
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internal state, but we are not free to choose a quantization direction. That is because the atom's dipole is induced 
by the vacuum fluctuations, and is free to point in any direction. In that light, choosing a particular direction seems 
invalid. Due to the form of the dipole - EM polarization function couplings, the induced atomic dipoles in different 
directions do not interfere, and a set of excited states (and thus different quantization directions) can be summed over. 
Such a set of independent excited states will form a resolution of unity and thus give a factor of unity contribution. 
The above transition amplitude can then be generalized to reflect the induced dipole: 

t+T 



(X /; t + T|exp[-~y ffds]|X i; t) 



DXexp 
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k 
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with p 2 = (g|p 2 |g) (the ground state expectation value of p 2 ). 



APPENDIX B: MOMENTUM COMPUTATION 



Putting in the spatial mode functions of Eq. © into the above gives the semi-classical transition amplitude in the 
presence of a conducting wall. 
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With the inclusion of the conducting boundary spatial mode functions the sums over momentum space are now 
over the positive half space. Despite it's complicated appearance, the transition amplitude above is in a useful form 
for computing the evolution of the momentum expectation value. The key point is that the transition amplitude 
of Eq. IBlfl is the product of several exponentials of exponentials, and contains only c-numbers. Therefore, each 
exponential can be expanded out into a series, the summands of all the series collected together, and the necessary 
integrations performed on the collected summand before redistributing the summand and resuming each exponential. 
That is, the individual exponentials in Eq. IjBlfl can be expanded in terms such as: 
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The resulting collected summand is 
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The momentum expectation value is then, 
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The momentum expectation value, the normalization factor, and other moments of the momentum operator can be 
computed with the generating function: 
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from which: 
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The factor A({n, n'}) is the summation measure. The initial wavefunction is taken to be a Gaussian centered at 
(R, Po) with the standard deviations (a,l/a). This choice allows the possibility that the atom is slowly moving 
toward the wall. Slowly, in this case, means adiabatically such that the external motion is much slower than internal 
time scales. 
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In the limits M — > oo and a — > such that & — > V and cr 2 M — > oo the generating function is: 



Z(J) = exp • 



/•t+T 

I ^ 

Jt 



iA 2 >' 
L 3 



cos 2 9 ; 



ik z -J 



e 2ik z .(R+Vs) _ ^ 



g 2 p 2 ' t+T 



_^ /-1-t-T /-S -, 



-i(w k +wo)(s-r)+ik-V(s-r) _|_ c c 



(B12) 



L 3 



ds fdrE^e^ 



e -i(w k +o;o)(s-r)+ik z -(2R+V(s+r-2t))+ik||-V(s-r) , c c 



^2 + iJ • Po + 0(e 4 /M) + 0{a 2 ) 



Finally, in the limits M — ► oo and a — * the momentum expectation value is: 
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